Monte Carlo Investigation of Ising Nanotubes and 
Nanostrips 



Carlos Garcia, M a Felisa Martmez and Julio A. Gonzalo § 

Departamento de Fisica de Materiales C-IV, Universidad Autonoma de Madrid, 
Cantoblanco 28049 Madrid, Spain 

Abstract. Monte Carlo simulations of the magnetization temperature dependence 
in D x L nanotubes (periodic lateral boundary conditions) and nanostrips (free lateral 
boundary conditions) with D=8, 16, 32, 64<C L — ► 5000 have been performed. The 
apparent critical temperature was determinated using the Binder Cumulant method 
(crossing of data for D x L with data for D x 2L) and it was found to be Tc = for small 
D values (D < D*), as it might have been expected. No scaling of ML^^ V vs |e| L x l v 
with \e\ = |T T J cl was fulfilled. However with |e| = e~ 4 / T (T c = for D < D* ~ 6), 

instead of |e| = T ~J°^ , scaling did hold. For D > D* the Binder Cumulant method 
produced clear phase transitions at Tq < Tc(d = 2), and susceptibility peaks were also 
observed. (Note that Tc is a non-universal quantity). The effective critical exponents 
did show a crossover towards one dimensional behaviour. The evolution of the apparent 
effective dimensionality (as defined by the effective critical exponents) for nanotubes 
and nanostrips as a function of D can be interpreted in analogy with what happens in 
thin films of thickness Dei: the effective exponents go from j3 e tf. m (D = D*) = 
to Peff. m (D » D*) ~ 1/8 and from S~^ m {D = D*) = to S~^ m (D » D*) ~ 1/15 
but, eventually, for L —> oo, and for T extremely close to Tc(D), both /3(D) and 
5^ 1 (D) cross over to j3(d = 1) = 0, 5~ 1 (d = 1) = 0, respectively as expected. However 
our Monte Carlo investigation did produce the novel result that there is a critical 
thickness, D = D* = 6 associated with a certain apparent effective dimensionality, 
which separates nanotubes with D < D* (for which Tc(D) = 0) from nanotubes with 
D > D* , for which Tc{D) is finite and positive. 



PACS numbers: 75.40.Mg, 05.10.Ln, 64.60.Fr, 05.50. +q 



§ To whom correspondence should be addressed (julio.gonzalo@uam.es) 



Monte Carlo Investigation of Ising Nanotubes and Nanostrips 



2 



1. Introduction 

Phase transitions in Ising systems with thin film geometry have been studied by many 
authors using series expansion, renormalization group and Monte Carlo calculations pQ- 
[H . Attent ion was focused on the evolution of the effective critical exponents at the 
vicinity of the critical temperature Tc, with the aim of determining the universality 
class of the thin films. The thickness dependence of the film critical temperature was 
analyzed by means of numerical finite size Monte Carlo simulations PU EI] , using free 
and periodic boundary conditions at the surfaces of the plates. It was found [T^HUI] that, 
for periodic boundary conditions, the Tc (D) scaling prediction for systems intermediate 
between d = 2 and d = 3, given by 



The scaling of the data and the crossover observed in effective critical exponents, 
which is due to the fact that the correlation length is sensing a finite thickness (D <C L) 
as the temperature goes towards Tc, demonstrated unambiguously that thin films belong 
to the two dimensional (d = 2) universality class. 

In this study we examine Ising systems with geometries intermediate between 
one dimensional and two dimensional, with the same purpose in principle, as in the 
investigation of thin films. Short-ranged isotropic Ising systems with a locally square 
lattices were investigated. In the case of d — 1-like systems it is clear that the 
determination of Tc, which is a non universal quantity, can be difficult as D approaches 
1, because the phase transition in a d = 1-infinite Ising system is attained only at 
Tc = and the exponents f3 and 5 _1 for the d = 1 Ising system are equal to zero in that 
limit. However this does not preclude the possibility that Tc > for beyond a given 
D* and effective critical exponents j3 e ff(D) > and 8~tAD) > appears which should 
eventually crossover towards (3 = and 5~ 1 = as T approaches Tc(D) very closely. 

Before presenting the results of our MC investigation let us anticipate what is 
expected to be analogous in the behaviour M(T) for strips of fixed width D and the 
behaviour M(T) for films of fixed thickness D, and what may be expected to be different. 
For strips we may expect that, at least for large enough widths D, there is a real phase 
transition with a Tc{L x D) approaching Tc(d = 2) = 2.269118... in the same way 
as for films of large enough Tc{L x L x D) was approaching Tc{d = 3) = 4.5115.... 
This effects the fact that even though strips of fixed finite D and length L — > oo belong 
certainly to the d=l universality class, the transition temperature is not a universal 
quantity and there for Tc may be zero for a while with D increasing from D = 2, 3, 4..., 
but it is not necessarily so for large enough D values. On the contrary for very large D 
values, Tc(D) must be expected to approach Tc(d = 2). We can expect in both cases, 
whenever T c > 0, that effective maximun critical exponents Peff.m and 8~^ m can be 
determinated for given D values, but also that for large enough L and temperatures 




(1) 



holds, where z/^ =3 ~ 0.63 is the correlation length exponent for three dimensional 
Ising systems. 
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very close to Tc{D) the effective critical exponents crossover to their one dimensional 
respective values ((3 = and = 0) for any D, in analogy with what happens with 

films ini-im. 

The difference between strips and films, however, is that in the former Tc(d = 1)=0 
and there is nothing that precludes Tc{D) = in a certain range of D values, f.i. 
D < D*. In contrast Tc(d = 2) = 2.269118 is non-zero from the beginning and films 
with small thickness D must be expected to grow from this Tq value from the very 
beginning. In the other words, Tc(D) a non universal quantity may have a threshold 
D = D* beyond which and only beyond which Tq{D) becomes positive and growing 
with D. 

2. Results 

In this work numerical finite size simulations of phase transitions in D x L Ising 
systems {D <C L) with free and periodic boundary conditions have been performed to 
investigate the width dependence of the critical temperature if any of strips and tubes. 
Metropolis Jl] and Wolff single cluster methods ^3] were used to perform Monte Carlo 
investigations of the temperature dependent magnetization M(T) in systems of width 
D=8, 16, 32, 64 and lengths 100 < L < 5000. Periodic boundary conditions were always 
used along the L direction to delete secondary boundary effects, while the conditions 
along the width (perpendicular to the L direction) of the system were modified from 
free to periodic in order to analyze boundary effects. Our analysis however concentrates 
on results obtained using periodic boundary conditions. 

The parameters of the algorithm (thermalization time, relaxation time and number 
of states for the energy count) were increased until finding that the main results did not 
change appreciably as these values were further increased and as it was confirmed that 
good statistics were obtained. The calculations involved 140000 Monte Carlo steps per 
spin for each data value (each temperature or each field). It may be noted that changing 
from 50000 MC steps to 140000 MC steps no noticeable improvements were observed. 
To reduce the critical slowing down as much as possible at the vicinity of the critical 
point we used a Wolff single cluster algorithm. Initial conditions at a given temperature 
were taken from the equilibrium conditions at the previous temperature. A description 
of the random number generator used may be found in ^BJ • 

Taking into account that Tc may be positive beyond a certain D* as explained 
above, we have investigated (i) the critical temperature dependence as a function of the 
width D, by simulating D x L systems to locate Tc {D,L — > oo), which for D larger 
than some D = D* is expected to be non-zero, (ii) the scaling of the magnetization 
data for two or more different L values with a constant D, using \Tj\ e = e~ 4// ' T_Tc ' for 
D < D*and e = ^ T ~J C ^ for D > D*, in the corresponding scaling variable. 

The apparent critical temperature for each strip (free boundary conditions) or 
nanotube (periodic boundary conditions) of width D, i.e T c (D, L — > oo) was determined 
using the Binder Cumulant [IE] Ul for two or more different L values (See below). We 
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did find that the crossing temperature between the Ul curves for a constant width D was 
strongly dependent of L for small D values. The L dependence of the Binder Cumulant 
Ul crossings for tubes and strips of width D is much more pronounced that the same 
dependence for thin films with thickness D jTUj but, eventually, as shown in Figure ^ 
below, it becomes stabilized for large L. Apart from possible corrections to scaling, this 
dependence on L may be connected with the fact that the Binder Cumulant is defined 
for second order phase transitions with an exponent (3 7^ 0, while our systems, for small 
enough D's, are approaching the linear chain critical behaviour, for which (3d=\ = 0. The 
value of Ul becomes increasingly difficult to determine exactly for smaller D's because 
the crossing occurs at the magnetization saturation region, but the crossing becomes 
increasingly an ambiguous for D D*. 

We can check the Ul crossings obtained with DxL and D x 2L, to define a "crossing 
critical temperature" Tc(D,L). This "crossing critical temperature" seems to approach 
a fixed Tq (D) 7^ as L — > 00, for D > D*. Figure^ shows that the "crossing critical 
temperature" obtained with the DxL and D x 2L systems for D > 6 tends to a fixed 
temperature value Tq(D) as L — ► 00 in such a way that < Tc (D) < Tq [d = 2) with 
T c (d = 2) = T C (D = 00) = \ln{l + sf2) = 2.269185314 as given by Onsager [TJj. An 
inset in Figure ^ shows a triple crossing for D = 64. 

The evolution of Tc{D,L) for increasing width D can be fitted approximately 
by means of an exponential. We take Tc(D), the transition temperature, as the 
extrapolated value at L — > 00. Table [T] gives the T c (D) values for periodic and free 
boundary conditions using the method described above. Again, we may note that Ta{D) 
is well defined but always somewhat overestimated in this way. It may be noted that the 
transition temperature Tc (D) for the same D value is higher for systems with periodic 
boundary conditions than for systems with free boundary conditions, i.e the nanotubes 
require more thermal energy than the strips to under go the phase transition, because 
their spins are more strongly bound due to the more stringent boundary conditions. 

It is true that the Binder Cumulant method, in particular the Binder Cumulant 
method crossing M(T) data for two or more finite size lattices, f.i. DxL, Dx 2L, Dx4L 
... is a phenomenological method but it is well grounded in general scaling arguments 
and, used with appropriate care, produces very reliable results. Figure [T] shows that it is 
necessary to go to relatively large nanotube and nanostrip lengths in order to stabilize 
the crossing temperature, but, of course, T C (D) begins to become very well defined as 
D grows beyond D ^> D*. 

As shown below to scale data properly for D < 6 it is necessary to use Tc (D) = 0, 
corresponding with Binder Cumulant but crossing at T = 0. We might say that D ~ 6 
corresponds to a certain intermediate effective dimensionality 1 < d < 2, as discussed 
below in more detail. 

Figure El gives (l — Tc(I^2) ) vs ^ ^ or P er i°dic (full circles) and free (open circles) 
boundary conditions in a log-log representation. The difference between both sets of 
points shows the effects of the boundary conditions on the critical temperature. An inset 
in Figure El shows a susceptibility peak for D = 64. Note that the scaling prediction for 
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systems of dimensions D x L is 

(l-^m oc D~ x ; A = — !— ; 1/^2 = 1 (2) 

and it does only hold for widths larger than D* ~ 6. This is in contrast to the thin 
films case [TU^ i n which the corresponding scaling prediction Equation (1) was fulfilled 
at any thickness with periodic boundary conditions. Thus, this plot serves to identity 
empirically D* = 6, below which all strips are effectively equivalent to linear chains for 
most practical purposes. 

Figure El shows a fit of Tc(D) 



T C (D) = T c (d = 2) [l - e~ mVD - D *J (3) 

to data for D = 16, 32, 64 (see Tabled) obtained with D* ta 6, m ta 0.353 ± 0.011. 
These T c values are used later for scaling purposes. 

We may note that the Binder Cumulant crossings are always close to 2/3. This 
means that the values found are almost equal to the d = 1 Binder Cumulant value. 
In all cases the values found are always far from the d = 2 Binder Cumulant value 
[20 . This fact implies that the universality class of the nanotubes (and nanostrips) 
must be d = 1. This fact does not preclude the existence of effective critical exponents. 
In order to investigate the scaling behaviour of nanotubes with D from D = 2 to 
D = 16 > D* = 6, we can use the following procedure: (a) for D < 6 we may take 
Tc = 0, and note that the appropiate scaling variable for e is e = e _4//T (identical 
to the case of the one dimensional chain). We note that (3d=x — and 1/Vd=i = 2 
did not scale properly for all D values. Then we did tray different values for f3 e ff 
and 1/ueff. /3 e ff = was correct for all D < D* and l/u e ff(D) was found to be 
l/z/ e// (2) ~ 0.897 ± 0.005, l/i/(4) ~ 0.371 ± 0.003. (b) On the other hand, for 
D > D* Tc(D) > 0, as given by Equation (3), and scaling was properly obtained 
with e = lT Tc(D) )l usin § &//(6) ~ 0.034 ± 0.001 and l/f e //(6) - 0.176 ± 0.005, 
ef f(8) ~ 0.0425 ± 0.001 and l/i/ e// (8) ~ 0.203 ± 0.003,and finally ,3(16) ~ 0.040 
and 1/V(16) ~ 0.247. This is illustrated in Figures 5-9. It may be noted that for D=32 
and D=64 the nanotubes (or nanostrips) begin to show that D is not properly D L 
with our data, L = 500, 1000, and the transition begins to appear less well defined, 
specially for the lower L value. 

Using the one dimensional chain scaling e = e~ 4 ^ T instead of e = y T ^ c j for 
D < D* , we can check that the exponent 7 e // evolves from j e ff.d=i — 1/2 towards 
r y e ff(D) > 1 as D increases. For D > D* , j e ff(D) decreases towards 7 e //.d=2 = 7/4. 

We note once more that for L — > oo and at T approaching very closely Tc(D) the 
effective exponents cross-over to the d = 1 exponents for any finite D <C L. 



3. Concluding remarks 



In conclusion, our Monte Carlo simulations show that: (a) Tc (D), which is a non- 
universal quantity, is Tc (D) = for D < D* = 6, and Tc (D) > 0,. It is well 
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approximated by Equation (3), being always smaller than Tq (d = 2) = 2.269185314, 
for D > D* = 6; (b) the universality class of nanostrips and nanotubes is one- 
dimensional but their phase transitions behaviour shows a qualitative change at the 
critical dimensionality D* = 6. Strips and tubes with very large D comparable to L 
are characterized by effective critical exponents approaching those of two-dimensional 
systems. 

It is clear that strips and nanotubes with D finite must be considered always one- 
dimensional systems. However our data show that regarding the transitions temperature 
T C (D), there is a definite change of behaviour at D < D* ks 6. A closer analysis of the 
dependence on D of the effective critical exponents is left for further investigation. 

It should be very profitable to investigate the temperature dependence of the 
magnetization by means of exact transfer matrix calculations, possible in principle for 
D values up to D > D* ~ 6. For D = 2 such calculations confirm unambiguously that 
the universality class is the same as for the linear chain (D — 1) with Tq{D = 2) = 
but it would be interesting to see what happens for increasingly higher values. Valuable 
information could be obtained in this way on the characteristics of the phase transition 
and perhaps on the functional form of the scaling function. This is left for further work. 

We hope to have shown here that Monte Carlo methods properly used, limited as 
they are, can provide very useful information on the phase transition features of low 
dimensionality Ising systems such as nanotubes and nanostrips. 

Acknowledgments 

We gratefully acknowledge financial support from DGICyT for grant BFM2000-0032, 
and a Universidad Autonoma de Madrid fellowship (C. Garcia). We want to thank C. 
Arago for help and advice. 



Monte Carlo Investigation of Ising Nanotubes and Nanostrips 



Table 1. Critical temperature (Binder Cumulant) for nanostrips and nanotubes. 
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Figure 1. Crossing critical temperature for D x L and D x 2L, as a function of 2L 
for nanotubes and nanostrips of width D=8, 16,32,64. The extrapolated temperature 
value at L — > oo is defined as Tc{D). The full and open circles correspond to 
nanotubes (periodic boundary conditions) and to nanostrips (free boundary conditions) 
respectively. The inset illustrates crossing of two Binder Cumulants L x D, 2L x D, 
4L x D. For L — > oo we could expect that the crossing take place at close temperatures. 



Monte Carlo Investigation of Ising Nanotubes and Nanostrips 

\D=6 



-a 

O 



1 F 



0.1 - 



0.01 



'0 



• Periodic 
Free 



D = 32 and L = K00O 



J V- 




2.25 2.50 



Slope ^-1.04 ± 0.05 



_l l l l_ 



T c (ld) 



10 



100 



D 



Figure 2. Width dependence of the critical temperature for nanotubes (periodic 
boundary conditions) and for strips (free boundary conditions) . The scaling prediction 
begins to hold for widths larger then D* ~ 6. The inset illustrates the fact that Tc > 
is accompanied by a pronounced peak in the susceptibility as a function of temperature. 
Note that the peak maximum appears at a temperature some what higher than the 
Binder Cumulant crossing temperature. 
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Figure 3. Transition temperature for nanotubes Tc(D) as a function of width D. 
Data from Binder Cumulants from Table I are well fitted by the equation show in the 
graph. 
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Figure 4. Effective critical exponent (3 e ff as a function of log \Tc(D) — T\ for various 
thickness. 
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Figure 5. Scaling plot of the magnetization for nanotubes of D=2 and L =500,1000 
using e = e~ 4 / T as scaling variable with Tc ~ 0. 
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Figure 6. Scaling plot of the magnetization for nanotubes of D=A and L =500,1000 
using e = e _4 / T as scaling variable with Tc ~ 0. 
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Figure 7. Scaling plot of the magnetization for nanotubes of D=6 and L =500,1000 
using e = ( ^Tc^ ) as scann g variable with r c «0. 
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Figure 8. Scaling plot of the magnetization for nanotubes of D=8 and £=500,1000 
using e = ( ^tI^ ) as scaling variable with Tc ~ 0.891. 
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Figure 9. Scaling plot of the magnetization for nanotubes of D=16 and £=500,1000 
using e — ( ^ T l r Tc ^ ) as scaling variable with T c w 1.525. 



